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In this paper we study the simplest massive 1+1 dimensional integrable quantum field the¬ 
ory which can be described as a perturbation of a non-unitary minimal conformal field theory: 
the Lee-Yang model. We are particularly interested in the features of the bi-partite entangle¬ 
ment entropy for this model and on building blocks thereof, namely twist field form factors. 
Non-unitarity selects out a new type of twist field as the operator whose two-point function (ap¬ 
propriately normalized) yields the entanglement entropy. We compute this two-point function 
both from a form factor expansion and by means of perturbed conformal field theory. We find 
good agreement with CFT predictions put forward in a recent work involving the present au¬ 
thors. In particular, our results are consistent with a scaling of the entanglement entropy given 
by ^ log I where Ceg is the effective central charge of the theory (a positive number related to 
the central charge) and I is the size of the region. Furthermore the form factor expansion of twist 
fields allows us to explore the large region limit of the entanglement entropy and find the next- 
to-leading order correction to saturation. We find that this correction is very different from its 
counterpart in unitary models. Whereas in the latter case, it had a form depending only on few 
parameters of the model (the particle spectrum), it appears to be much more model-dependent 
for non-unitary models. 
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1 Introduction 


Entanglement is a fundamental property of quantum systems which relates to the outcomes 
of local measurements: performing a local measurement may affect the outcome of local mea¬ 
surements far away. This property represents the single main difference between quantum and 
classical systems. Technological advances have taken entanglement from a strange quantum 
phenomenon to a valuable resource at the heart of various fields of research such as quantum 
computation and quantum cryptography. There has also been great interest in developing ef¬ 
ficient (theoretical) measures of entanglement, not just in view of the applications above but 
also as a means to extract valuable information about emergent properties of quantum states of 
extended systems. One such measure for many-body quantum systems is the bi-partite entan¬ 
glement entropy (EE) [T], which we will consider here. Other measures of entanglement exist, 
see e.g. [n El la n E], which occur in the context of quantum computing, for instance. In its 
most general understanding, the EE is a measure of the amount of quantum entanglement, in 
a pure quantum state, between the degrees of freedom associated to two sets of independent 
observables whose union is complete on the Hilbert space. In the present paper, the two sets 
of observables correspond to the local observables in two complementary connected regions, A 
and A, of a 1-|-1-dimensional (1 space -|- 1 time dimension) extended quantum model, and we 
will consider cases where the quantum state is the ground state of a non-unitary, near-critical 
model. 

Prominent examples of extended one-dimensional quantum systems are quantum spin chains. 
Their entanglement has been extensively studied in the literature laEiEiEiiiniiiiiiiaiisKii]. 
These examples however all refer to unitary quantum spin chains. Interesting examples of non- 
unitary spin chain systems exist, for instance the famous quantum group invariant integrable 
XXZ spin chain, with generically non-Hermitian boundary terms; in the thermodynamic limit it 
has critical points associated with the minimal models of conformal field theory (GET), including 
the non-unitary series [laESKITKls]. Another example is provided by the Hamiltonian studied 
by von Gehlen in |19[ EO] : the Ising model in the presence of a longitudinal imaginary magnetic 
field. This Hamiltonian has a critical line (in the phase space of its two couplings) which has been 
identified with the Lee-Yang non-unitary minimal model of GET, with central charge c = —22/5 
[SUES]. In all these examples, the local, extended Hamiltonians are non-Hermitian, yet have 
real and bounded energy spectra. Their critical points are described by GET models containing 
non-unitary representations of the Virasoro algebra with real weights, and whose ground states 
are not the conformal vacua, but negative-weight modules. 

Non-Hermitian Hamiltonians with real spectra are the subject of much current research 
especially in connection with PT-symmetry or pseudo/quasi Hermiticity |23l El] (see |25l E6l EZj 
for reviews and [28] for the interplay with integrability). For instance the critical line of von 
Gehlen’s system HaED] described above, can be related to PT-symmetry breaking in that it 
separates the phase space into two regions, one where only real eigenvalues occur, and another 
where pairs of complex conjugated eigenvalues arise [29]. Experimental studies and theoretical 
descriptions of new physical phenomena connected to non-Hermitian Hamiltonians have recently 
emerged, including optical effects [SniEIlES], transitions from ballistic to diffusive transport 
|33| . and dynamical phase transitions |34L I35| . Non-Hermitian quantum mechanics is also used 
in the description of non-equilibrium systems [36], quantum Hall transitions m, and quantum 
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annealing [38]. 

At quantum critical points, the scaling limit of the EE has been widely studied within 
unitary models of GET [39l SOI [T] E] HU H2] . In particular, the combination of a geometric 
description, Riemann uniformization techniques and standard expressions for GET partition 
functions is very fruitful. Recently [H], this was generalized to non-unitary GET, where a 
general formula was obtained using such techniques. Near critical points, the scaling limit is 
instead described by massive quantum held theory (QET), and geometric techniques relying 
on conformal mappings break down. As was found in [44114^ 146] . the most powerful way of 
studying the EE in unitary models of QET is using an approach based on local branch-point 
twist fields. However, the question of the EE in non-unitary near-critical models is much more 
delicate, and standard arguments give little indications as to how to modify the held-theoretical 
approach. Importantly, the rigorous derivation presented in |l3j provided a precise local-held 
description of the EE involving composite helds in the branch-point twist family, thus opening 
the door to its study in non-unitary QET. In the present paper, using techniques of integrable 
QET, we will study the scaling limit of the EE in the near-critical region of von Gehlen’s model, 
described by the Lee-Yang QET model. 

This paper is organized as follows. In sectionwe recall the main dehnitions and techniques, 
and provide a summary of our main results. In section we introduce the Lee-Yang model and 
some general results on the form factor expansion of correlation functions, their logarithms and 
expectation values of local helds. In section]^ we review the twist held form factor equations and 
present solutions for the branch-point twist helds helds T and : T(p : in the Lee-Yang model. 
In section we test our form factor solutions by performing a form factor expansion of the 
functions log ^(: T4> :)“^(: T4> : (r) : T(p : (0))^ and log (^(T’)“^(T(r)T(0))^ and recovering the 
behaviours log(mr) and —Axq- log(mr) for some constants x-p, which we compare to 

GET predictions. In section]^ we compare a form factor computation of the two-point functions 
above with a computation in zeroth order conformal perturbation theory. As a byproduct, we 
hnd general formulae for some of the GET structure constants entering the OPEs of T and T 
and of : T(p ■ with :'T(p'-- In section]^ we present numerical results for the Renyi entropy near 
criticality and a detailed computation of the hrst three leading corrections to saturation of the 
EE. We hnd that the next-to-leading order correction to saturation is non-universal. In section 
[^we present our conclusions and outlook. In appendixwe explain how the normalization and 
conformal dimension of the held : T(p ■ are hxed by GET. In appendix]^ we present a detailed 
analysis of the one-particle form factor contribution to the two-point functions of T and :Tcp 
A large n expansion of this function demonstrates that it provides a very substantial contribution 
to the power law behaviour of the two-point functions at short distances. In appendix we 
present a computation of the three particle form factor of Lee-Yang twist helds. In appendix 
l^we perform a computation of some of the structure constants entering the OPE of helds T 
and T and of helds : T(p '■ and : Tcp : in GET. In appendix we present a computation of 
the numerical coefficient of the next-to-leading order correction to saturation of the EE in the 
Lee-Yang model. 
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2 General aspects and summary of main results 

In order to provide a formal definition of the EE, the Hilbert space of an extended quantum sys¬ 
tem, such as a spin chain, is decomposed into a tensor product of local Hilbert spaces associated 
to its sites. Grouping together sites associated to the regions A and A, this gives: 

n = A®A. ( 1 ) 

The EE in a state \'ip) is the von Neumann entropy of the reduced density matrix pA associated 
to A: 

Sa =-Ti'APA'^ogpA , M = Tr_4|V^)('0|. (2) 

Another frequently used measure of entanglement is the Renyi entropy, 

oA) _ logTV^pl 

which specializes to the von Neumann entropy at n = 1, 

lini = - lini ^TvaPa = Sa- (4) 

We will study the ground state entanglement entropy in the scaling limit of infinite-length 
quantum chains. The scaling limit gives the universal part of the quantum chain behaviour 
near quantum critical points, described by 1-|-1-dimensional QFT. It is obtained by approaching 
the critical point while letting the length i of the region A go to infinity in a fixed proportion 
with the correlation length | (measured in number of lattice sites). If ^ = oo from the start, 
the system is exactly at its critical point, and the scaling limit is described by GET. In this 
case the entanglement entropy of unitary critical systems, as a function of i, is divergent in 
a way which was first understood in [391 l4Uj . numerically confirmed in mu and generalized 
and reinterpreted in ED ig. The divergency is logarithmic with a proportionality constant 
depending on the central charge c of the GET, 

= + SA{(.) = l\og-^+o{l) (GET), (5) 

and where e is a non-universal ultraviolet cut-off (proportional to the lattice spacing) which is 
chosen so as to encode all o(l) corrections. The formulae above are easily adapted to the case of 
an inhnite region ^ = oo near criticality ^ < oo, where i is simply replaced by ^ in @ ED 132]. 
In the full scaling limit, where (. and ^ are both large and in proportion to each other, there is 
a universal scaling function which interpolates between the two results, 

S^a\^) = log - + /W6 + 0(1) (QFT). (6) 

bn £ 

In this case the result is much less trivial and has been studied in unitary integrable [441 05] and 
non-integrable ES] models using massive QFT techniques. 
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We may ask how (if at all) the entanglement entropy is affected by non-unitarity. At criti¬ 
cality, it was shown in [33] that the entanglement entropy scales instead as 

^ log ^ -b o(l), Sa{() = ^ log ^ -b o(l) (non-unitary CFT), (7) 

where c^s ■= c — 24A, and A is the smallest (often negative in non-unitary models) scaling 
dimension of a primary held in the CFT. For the Lee-Yang model, for example, Ces = | as 
A = — This result is not entirely surprising as the work of Itzykson, Salem and Zuber m 
had previously shown that the effective central charge Cefr also replaces c in the expression of the 
ground state free energy found by Affleck [38| and Blote, Cardy and Nightingale [39]. However, 
the question of the entanglement entropy in non-unitary near-critical models is much more 
delicate. Importantly, the rigorous derivation of Q presented in |43j has lead to new insights 
into the computation of entanglement entropy in non-unitary theories and its held theoretical 
interpretation, opening the door to its study away from criticality in QFT. 

It is known since some time [Ml SQi im sa that the bi-partite entanglement entropy in the 
scaling limit can be re-written in terms of more geometric quantities, using a method known 
as the “replica trick”. The essence of the method is to “replace” the original QFT model by a 
new model consisting of n copies (replicas) of the original one. These are used to represent 
when n is an integer, and then to evaluate Tyap'a- The quantities Sa and 5^^ for general n are 
then obtained by “analytic continuation” in n. The matrix multiplications in and the trace 
operation give rise to the condition that the copies be connected cyclically through a hnite cut 
on the region A. As a consequence, this trace is proportional to the partition function Z„(xi, X 2 ) 
of the original (euclidean) QFT model on a Riemann surface A4n,xi,x2 with two branch points, 
at the points xi and X 2 in and n sheets cyclically connected. The positions xi and X 2 of the 
branch points are dimensionful positions in the QFT model corresponding to the end-points of 
the region A in the scaling limit. This gives: 


5a (r) = — lim 


d Zn{xi,X2) 


n->i dn 


Zf 


( 8 ) 


Here, r := |xi — X 2 I is the euclidean distance between xi and X 2 - The above concepts hold, 
in principle, for any QFT model, unitary or not. In CFT, one may evaluate this by using 
the uniformization theorem: the Riemann surface Ain,xi,x 2 can be conformally mapped to the 
Riemann sphere with two punctures (or the cylinder) by using the map g reproduced in appendix 

El 

In the EE context, it was first noticed in |41ll42j that the ratio of partition functions above can 
be reinterpreted as correlation functions of certain fields, which were not otherwise specified, 
in unitary CFT. This idea was then generalized to unitary massive theories in [33] and the 
fields where identihed as branch-point twist fields T(xi), T{x 2 ) characterized by their non-trivial 
exchange relations with other fields of the n-copy theory. These twist fields are defined only in 
the replica model (e.g. they become the identity field when n = 1), and are primary fields arising 
from the extra permutation symmetry present in the replica theory; they are associated to the 
Z„ symmetry generators j 1 —)• j -b 1 mod n and j 1 —>■ j — 1 mod n respectively. In CFT, such twist 
fields and their relation to partition functions on Riemann surfaces were in fact studied much 
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before their use in the computation of the EE was emphasized, see for instance m- In terms 
of these fields, the replica partition function is given by 

= Zr.e^^r(^r{x^)f{x2)), (9) 

where {T{xi)T{x 2 )) is a two-point function in the ground state of the replica theory. The branch¬ 
point twist fields are chosen so as to have the GET normalisation (e.g. the leading term in their 
OPE has coefficient 1). The constant Zn, with Zi = 1, is an n-dependent non-universal constant, 
e is a short-distance cut-off which is scaled in such a way that dZn/dn = 0 at n = 1, and, finally, 
Aj- is the conformal dimension of the counter parts of the fields T, T in the underlying re-copy 
conformal field theory, 

which can be obtained by GET arguments |5n[ W2\. ITS] . It is easy to show that the formula Q 
when inserted in Q indeed reproduces ([^ for GET. 

The derivation above assumes unitarity of the theories under consideration. In such case Aj- 
is by construction the lowest conformal dimension of any field in the replica theory which has the 
twist property. The GET derivation of ([^ has been generalized to the non-unitary case in [l3| 
leading to the expressions Q . In this work it was also observed that the EE could be computed 
from a representation of the Z^-orbifold partition function of the theory via correlation functions 
involving certain twist fields of the re-copy replica theory. This representation requires new 
twist fields : Tcj) : and : T^(p obtained from the primary twist fields T and T as leading 
descendants in the product with the lowest-dimension field cj) (of conformal dimension A). More 
precisely: 

n 

: T(j) ■ {y) = lim |x - ^T{y)(j)j{x), (11) 

i=i 

and similarly for :'T<p:. These composite fields were first introduced in EH and further studied 
in [52|. The constant ensures conformal normalization, namely 

(: Tct> : (xi) : fcj) : (xa)) = |xi - (12) 

in GET, where 

(13) 

are the conformal dimensions oi : T4> '■ and : X4> In the context of the study of the EE they 
were hrst obtained in m- However, as for many other quantities in this context, they had 
emerged previously in the study of orbifold GET, see e.g. [MlEaiSS]- 

A detailed derivation of both the normalization constant and the power law in 0 is given in 
appendix [Aj The dimension A.-y^. arises naturally in computations of the entanglement entropy 
in non-unitary GET, and, as noticed in |l3|, suggests that for such theories, the partition function 
of the re-copy theory may be written instead as: 




Zn{xi,X2) ^ {■ T(p : (xj) : t'(/> : (X2)) 

" (<()(X1)</-(X2))- 


(14) 
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Figure 1: An artist’s impression of the Riemann surface Mn,xi,x 2 for n = 3 with held insertions 
(p at the branch points xi, X 2 - 


where again Zn is such that it and its derivative at n = 1 are 1 and e is a short distance 
cut-off. Compared to Q the expression (14) involves not only a different twist held but also 
a normalization given by {(j){xi)(j){x 2 ))'^■ For CFT it is easy to interpret this normalization as 
simply the norm of the ground state which in radial quantization is created by the action of the 
held (p on the conformal vacuum. As for the unitary case, it is easy to show that the formula (14) 
when inserted in (Q reproduces ([^ for CFT. It is natural to assume that the same expression 
will hold beyond criticality. This paper is a hrst step towards putting this assumption to the 
test beyond criticality. 


2.1 Summary of main results 

From the formulae above it is clear that a study of the EE in massive QFT is in principle 
only possible by studying correlation functions of twist helds. This approach has been pursued 
successfully in several works |44l [56l 145] where the ratio of partition functions (|^ at large 
distances r = |xi — X 2 I (the infrared (IR) region) has been studied for unitary 1-|-1-dimensional 
integrable QFTs. Integrability means that in these models there is no particle production in 
any scattering process and that the scattering (S') matrix factorizes into products of two-particle 
S-matrices which can be calculated exactly (for reviews see e.g. [571 EH ESI ESI El])- Although 
most of the integrable theories studied in this framework are unitary, well-known examples of 
non-unitary integrable QFTs exist. Best known among those examples is the Lee-Yang model 
whose exact S-matrix was first given in |62j . 

Taking the known S-matrix of an integrable model as input it is possible to compute the 
matrix elements of local operators (also called form factors). This is done by solving a set of 
consistency equations |631 [64] , also known as the form factor bootstrap program for integrable 
QFTs. It is this particular feature of integrable models which makes them interesting, as it 
means there is a systematic, non-perturbative way of computing multi-point functions of local 
helds. These computations are by no means easy, but often provide good numerical accuracy. 
In |44j . the form factor program was generalised to branch point twist helds leading to the 
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evaluation of Q for various unitary models. In this paper we will pursue this program for the 
Lee-Yang model employing the formula (14). Our main results can be summarized as follows: 


1) Form factor program for twist fields: We have found that the twist field form factor 
equations together with the requirement of form factor clustering are sufficient conditions 
to entirely fix all form factor solutions for any particle numbers. In particular, these 
constraints immediately give rise to two form factor families, naturally identifiable with 
the fields T and : Tf) We have carried out a zeroth order perturbed CFT computation of 
the twist field two-point function for several values of n and compared this to a truncated 
form factor expansion of the same correlator. The former is expected to be accurate at 
short distances, the latter at large distances. Nevertheless, the agreement is relatively 
good, thus confirming the validity of the form factors found. 

2) Saturation of the EE at large subsystem size: Let us absorb all non-universal o(l) constants 
of the short-distance behaviour of the Renyi EE into a short-distance cutoff e„. Subtracting 
this non-universal contribution, the EE at large distances then saturates to a universal 
constant which can be calculated using QFT. More precisely, we find 


~ + 1) + Un + o(l) {mr oo) 


6n 

Ces{n + 1 ) 

6n 


log(r/en) -h o(l) 


(mr —)• 0) 


where the universal saturation Un is given by 

1 

U — ^ 

^ n. — 


1 - n K,rcl>: 

The constants Ko are fundamental properties of QFT fields O, defined by 

lirnmr^o (mr)^*o(C)(r)O(0)) 


Ko = 


limn 


(O(r)O(O)) 


(15) 


(16) 


(17) 


where xq is the unique exponent making the limit finite and nonzero, and where O is the 
“conjugate” under internal symmetries {(f) = f and : T4> ■ =■ 'T :). In the unitary case, 
Xq is the conformal dimension of O, and Kq = {0)~^ under the CFT normalization 

of O. In the non-unitary case, both of these statements are modified. In particular, in the 
Lee-Yang model, these constants can be expressed in terms of massive QFT and of CFT 
data as 


K^ 


= m 


, 2 A 



(<)>)’ 


K,r4>: = 


(: Tf :)2 


(18) 


where the vacuum expectation values are under CFT normalization, and and 

are structure constants of conformal OPEs. Constants Kq can also be expressed solely in 

terms of form factors of the massive model, as in (32). 
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3) Leading order correction to saturation: For unitary theories, one of the most interesting 
results |44l [56] has been the identification of a universal leading order correction to the 
large-distance (large-r) saturation of the entropy of all unitary integrable theories. This 
exponentially decaying correction, of order o{e~‘^^^), has a higher degree of universality 
than usual QFT quantities, as it only depends on the particle spectrum of the model. In 
|46j it was shown that, even more strikingly, this feature holds beyond integrability. This, 
however, seems to be broken in non-unitary models. For the Lee-Yang model, we found 

SA{r) ^ log(me) + U - aKo{mr) + (19) 

15 

where a = —0.0769782... is a constant that is (a priori) model-dependent. This suggests 
that we may use this feature of the EE as a means to identify non-unitary critical points. 
Indeed, given a spin chain model whose critical point is not known, a study of entangle¬ 
ment at criticality will reveal the value of c^s- However, this does not say if Cefr = c or 
not. Considering large size corrections away from criticality will reveal different types of 
exponential decay depending on whether or not the theory is unitary. 


3 ^-matrix and form factors in the Lee-Yang model 

3.1 S'-matrix 


The Lee-Yang model is one of the simplest 1-1-1 dimensional integrable QETs. Erom the GET 
point of view, it may be regarded as a perturbation of the non-unitary minimal model associated 
with central charge c = — The primary operator content of the theory is very simple, 
consisting of the identity and a scalar field </> of conformal dimension A = — g. Perturbing this 
CFT by the scalar field we obtain the massive Lee-Yang model. This theory has a single particle 
spectrum. The scattering amplitude corresponding to the scattering of two particles of the same 
type was found by Cardy and Mussardo [^2] and can be written as 


S{9) 


tanh 5 (0 + 
tanh \ 


It has a pole in the physical sheet at 0 = ^ corresponding to the formation of a bound state, 
which in this case is the same fundamental particle of the theory. We note that the non- 
unitarity is manifested by the fact that the associated residue has the wrong sign. Nevertheless, 
the corresponding integrable massive model is well defined. The n-copy model, where the Z„ 
twist fields live, possesses n particle species /r = 1 ,... ,n, and a two-particle scattering matrix 
given by = 5 ( 6 »)'^^‘i''^ 2 . 


3.2 Form factor expansions of two-point functions 

In this paper we study the correlators (Y(r)T(O)) and, especially, {: T4> : (r) : T(p : (0)) and 
((^(r)())(0)), as well as the associated entanglement entropy obtained via (14). It is well known 
that two-point functions of local operators in QFT can be expressed as infinite sums involving 




matrix elements of these operators. The matrix elements of relevance, also known as form 
factors, are defined as 

Ff'^-^'=(01, ...,ek):= (O|O(O)|0i,..., , (20) 

for a local field O. Here |0) represents the vacuum state and |0i,..., the physical 

“in” asymptotic states of massive QFT. They carry indices m, which are quantum numbers 
characterizing the various particle species, and depend on the real parameters 9i, which are 
called rapidities. The energy and momentum of a particle of mass rm are expressed in terms 
of its rapidity 9i as mi cosh 9i and mjsinh^j, respectively. In terms of form factors, two-point 
correlation functions (in unitary models) may be expanded as 


(OMO'(O)) 


oo 


E 


1 


E 





{01, ■■■,6k) 


k 

2 —rm cosh^j 

e 


( 21 ) 


As mentioned, the Lee-Yang model is non-unitary. As was noted in |65[ 166] . non-unitarity 
affects the form factor expansion. A consequence of this is that many fields appear to be non- 
Hermitian under the Hilbert structure of asymptotic states. In the Lee-Yang model, an exact cal¬ 
culation of the form factors of the field cj) shows that ((0|())(0)|6li,..., 0^)“)* ^ • • •, 0A:|'/>(O)|O), 

where the right-hand side can be obtained by crossing symmetry. However, it turns out that 
the relation is surprisingly simple: 

((</.)-i (O|^(O)|01,..., = (-1)" (0)-i “(01,..., 0fc|0(O)|O). (22) 


As a consequence, the form factor expansion of the two-point function of the field cj), normalized 
by the square of the VEV, is a modification of (21) where sign factors (—1)^ are included for the 
terms involving the /c-particle form factors. This gives rise, in the single-copy Lee-Yang model, 
to: 


((/.(r)(/.(0)) 

(</,)2 


E 

k=0 


(- 1 )^ 


k\ 





—rm cosh^j 

e 


(23) 


A natural way to understand this modification is through a discussion of the bound-state 
singularity occurring in the form factors. The additional (—1)^ guarantees that the bound- 
state residue of the analytic continuation of the /c-particle integrand, which, like that of the 
scattering matrix, has the wrong sign, is related to the k — 1-particle integrand in a way that 
would guarantee locality properties. As we will see below, form factors of twist fields : 'T4> '■, 
: T(j) '■, T and T are subject to similar bound-state residue equations as those of 4>. Hence, this 


interpretation suggests that a similar modification of (21) occurs for the form factor expansion 


oi {: T4> : (r) : T(p '■ (0) and of (T(r)T(O)). That is, in the n-copy model. 


(0(00(0)) 
{O? 



— 1 77'^IMI 


(0)-^F 


(01, ...,0k) 


—rm cosh^j 

3 0=1 


(24) 
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for both O = T, O = T and O =: T4> O =: Tcp Here we have used the fact that by 
symmetry under inversion of copies, (T) = (T) and {: T(p '■) = {'■ T(p :), and the form factor 
expansion includes sums over the copy numbers fij. That this is the correct expansion follows 


from an equation similar to (22) for twist fields, see subsection 4.4 


Finally, since the field p is no longer Hermitian, its VEV is no longer expected to be real. 
As was shown by Zamolodchikov, {p) is in fact purely imaginary - this can be explained by the 
fact that it occurs in the formal massive Lee-Yang action (written as a perturbation of the CFT 
action) with a purely imaginary coupling constant. A similar phenomenon makes the VEVs 
(T) and (: Tp :) not necessarily real. We will determine their phases (up to multiples of tt) 
by evaluating analytically the normalization of their leading short-distance power-law, and by 


observing numerically that the right-hand side of (24) is positive for all mr. 


3.3 Short-distance behaviour from form factors 


Form factor expansions (21) and (23) are naturally large-distance expansions, in that they 


converge very rapidly for large values of rm. However, in many cases we want to explore small 
values of rm. In such cases two-point functions generally develop power-law behaviours in rm 
and it is very difficult to extract the precise power from a form factor expansion such as those 
above. 

It was realized a long time ago [65] (see also ISZI for a nice derivation and application to 
various models and [68] for a generalization to boundary theories) that if one is interested in the 
short-distance behaviour of correlators then an expansion of the logarithm of two-point function 
is more appropriate: 





(Oi,- 


—rm cosh^j 






.(25) 


The functions 

recovered when exponentiating (|25|). This condition automatically implies for example that 


,dn) must of course be chosen so that the expansion (21) is 






TjO\iiifi2^3 {n n A 


(C>)■Vf'^'(^)l^ (26) 

- Hfl^'(0i)Hfl'‘"(02), (27) 

(28) 


In general the Hj. functions can be interpreted as the “connected parts” of the functions 
(they are “cumulants” with respect to the rapidities). These are such that, if the clustering 
decomposition holds for the F^s at large rapidities for all k, that is 


lim F, 

ei,...,6»fe->-oo 


,0|/ii fp a \ ^k (.I'l) ■ ■ ■ 1 ^k)r^ (^pfe-i-i, . . . , C/fc-|_^j 

k+e [^ 1 , ■ ■ ■ ,^k+e) — ; (29) 
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\^k,^ G N, then the HkS vanish at large rapidities for all k. 

Thanks to this vanishing, for mr <C 1 we now expect each summand in the sum over k in 
the expression above to be dominated by a leading term proportional to log(mr). The constant 
coefficient of this term, summed over all particle contributions, will then give the power which 
governs the short-distance behaviour of the two point function. Let us call this power —Axq. 
Then, carrying out one integral in (25) and expanding the result for small mr we find 

k 


1 ^ (- 1 )^ 

= — 2 ^ 


dvr 


k=o 


k\ 


s n, 

Ml V 1 \.j —2 _ 





( 0 , 02 ,• • • , 9k). 


(30) 


This was used for the field (p in | 66 ] and shown to agree well with conformal field theory results. 
In addition, the proportionality constant (0 of the power law behaviour at short distances. 




[mr 


0) (31) 

can also be extracted from the form factor expansion. It was shown in |67! that by considering 
the leading correction to the log(mr) term in (25) one may also find a form factor expansion for 
the constant Kq which is given by: 


iLo = exp I 


k=0 


k\ 


with cosh 0 j + Ij - (^X]j= 2 sinh 0 j^ 

Mascheroni constant. 





(0,02,-•• ,0fc)(ln|+7) 1 , (32) 


and where 7 = 0.5772157... is the Euler- 


4 Twist field form factors 

4.1 Form factor equations and minimal form factors 

The form factor equations for Z„ twist fields were derived in |44j . Details of the solutions 
procedure for two-particle form factors appeared there, and higher particle form factors of various 
models were computed in |45] and |69| . Interestingly, a very similar set of form factor equations 
had been derived much earlier m in a rather different context (e.g. the study of the response of 
an integrable QFT to a variation of the Unruh temperature). The details of the computation for 
the Lee-Yang model are very similar to those described in these works, with the only difference 
that the presence of the bound state pole in the 5-matrix imposes further conditions on the 
form factors. In particular, bound state poles are present in addition to kinematic poles. The 
form factor equations only encode locality properties of helds, hence they are unchanged for 
form factors of any held in a the same twist sector. In order to distinguish for form factors 
of T and : T 0 :, we will impose additional conditions, and verify the correctness of the solutions 
by numerical comparisons with CFT predictions. 
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In what follows we will consider /c-particle form factors 62 , ■ ■ ■, Ok) for a generic 

twist field O. We will later identify this field with T or : Tcj) '■ depending on various properties 
of the form factor solntions we obtain. 

The two-particle form factor mnst satisfy (in the two-particle case we nse the single argnment 
0 = 61 - O 2 ) 

+ 2mn), (33) 

and the kinematic residue equations 

Res Ffl^^(0 + i7r) = i(0), (34) 

0 = 0 

Res F^^^^{0 + m) = -i{0). (35) 

0 = 0 

Here and below we use jX = ^ — 1 mod n. Higher particle versions of these equations read 


Res 

00 = 00 

■■^''(0o-Fi7r,0o,0i.. 

■ , 0 k) 


(36) 

Res 

00 = 00 

■■^'=(0o-Fi7r,0o,0i.. 

■ , 0 k) 

i=l 

,.,0fc). (37) 


For this model there is the added difficnlty of having to solve also the bound state residue 
equation associated to the scattering process a + a ^ a where a is the Lee-Yang particle on copy 
a. This takes the form 


Res F' 




= 0 


n+1 


(0 + y, 0 - y, 01, • • ■, 0n-i) = iTFf (0, 0i,..., 0„_i) (38) 


where the so-called three-point coupling is fixed by 

^TT? 

= —i lim (0-)5'(0) = — 2 -V /3 (39) 

3 

and by choosing the negative imaginary direction: T = —For n = 1 this eqnation fixes 
the one particle form factor (which for spinless fields must be rapidity independent) through the 
equation 


Res l““(0 - 0 + —) = iFTf (40) 

0 = 0 ^ 


These equations imply that the two-particle form factor solution given in [H] mnst be gen¬ 
eralized to include the bound state pole. As discussed in [66] this may be done by defining a 
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minimal form factor. A minimal form factor is a solution of (33) which has no poles in the 


(extended) physical sheet 0 G [0, 27rn) except possibly for bound state poles, and which tends 
to unity as |0| —)• oo. It turns out that this particular Riemann-Hilbert problem has a unique 
solution, and this solution possesses bound state poles with nonzero residues: 


Fmm{0) = a{6,n)f{e,n), 
where a{9, n) encodes the bound state pole 


a(0, n) = 


cosh - — 1 

_ n _ 

cosh ^ - cos ^ 


(41) 


(42) 


and f{9) is given by the integral representation 
f{9, n) = exp 1 2 


°° sinh I sinh | 
3 6 


0 t sinh(nt) cosh | 


cosh t { n + 


i9 


TT 


(43) 


The latter function admits also a representation as an infinite product of gamma functions which 
was already given in |44] for the sinh-Gordon model (it suffices to take i? = 2/3 and to invert 
the formula). 


The expression (43) may be obtained as a solution to (33) using a similar integral representa¬ 


tion of the two-particle scattering amplitude. In the absence of bound state poles, the resulting 
f{9, n) would directly be the minimal two-particle form factor. In the present case, however, 
the function tends to 1 as |0| —)■ oo but has a simple pole at 0 = 0. The factor a(0,n) is the 
unique one that shifts this pole towards the position of the allowed bound-state singularity in 
the physical sheet, without affecting the large-|0| behaviour. 

Using the integral or Gamma-function representation, it may be shown that 


n 


sin" ik 


/(^7r,n) _ 

/(^,n)2 VSsin^sin^' 


(44) 


In order to compute higher particle form factors the following more general identities are im 
port ant 

cosh ® - cos 




3n 


.(0+f)i^min(0-f)= , 

o 3 cosh 5 — cos ^ 


Fmm(9 + Z7r)Tmin(0) = 


™l>|;Sinh(|; + g) 


(s - s) ( 


0 I 5i7r ) 
2n ' 6n J 


(45) 

(46) 


4.2 Twist field one- and two-particle form factors 

One expects that primary twist fields, with direct geometric meaning, would occur as solutions 


to (33) and to conditions of bound state and kinematical singularities with the additional re¬ 
quirement of convergence as |0| —)• oo. This additional requirement is expected to implement, in 
a path-integral picture, the least singular asymptotic condition possible at small distances near 
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the position of the field. With these conditions, the most general form the two-particle form 
factor can take is 


F?"\e) = 




FraUO) 


2nsinh (^) sinh (^) Fmin(*vr) 




(47) 


where the first term is of the form required to solve the kinematic residue equation (and of the 
same form as for other theories previously studied [H] ) and the second term is what is commonly 
termed a “kernel” solution of the kinematic residue equation (that is a solution without kinematic 
poles). 

In general k is an arbitrary constant, but it may be fixed by imposing the cluster decompo¬ 
sition property, namely 


lim F^'^^(0) = K : = 


6^00 


(o) ■ 


(48) 


where we have used the fact that limg_,.oo (9) = 1. Then, the one-particle form factor on 
copy a may be fixed by combining this with equation (40), which translates into the following 
quadratic equation for 


C.|1_ 1 tan s/(¥.»),p,, 

--rtani /{„,„) + 


TT 

3n 


,,27ri 


This leads to two possible solutions: 


Ff = -{0)T 


cos(l^) ±2sin^ (^) 
2’^sin(|^)/(^,n) ’ 


(49) 


(50) 


where we have used the identity (44). 

The presence of two solutions immediately suggests the existence of two different least- 
singular twist fields, by contrast to other models studied in the past. It is natural to conjecture 
that these are 7~ and : T4> and given this, it is a simple matter to identify their respective 
form factor solutions. Indeed, the former specializes to the identity at n = 1, and the latter, to 
(j)- We note that the solution with the negative sign specializes to 0 at n = 1, and that that with 
the positive sign specializes to the one-particle form factor of the field (p 


i2V2 


3i/4/(^,l) 


with Fq = {(p) = 


him 


24/iv/3 


and h = 0.09704845636... 


(51) 


found in [66l (note that the constant u(0) in |66] is u(0) = /(zvr, 1)^/^ = ^/(^,1) with our 
present notation, and that the coupling h was computed in m)- These properties suggest the 
identifications 




-r|i 


= -T- 


2cos(^) -1 




:¥</>: 11 


(T) " 2nsin (^)/(^,n)’ {■ Tp :) 2resin (^)/(^, n) ‘ 

The numerical results of sections 4, 5 and 6 provide further support for these identifications 


(52) 
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4.3 Higher particle form factors 


Let us now consider only form factors of the form , Xk) '■= (xi,..., Xk), that 

is form factors involving only one particle type. This is sufficient as form factors involving other 
particles may be obtained from these by using the twist field form factor equations |44j . 

The higher particle form factors may be obtained by making the ansatz 


F^{xi,...,Xk) = Qkixi,...,Xk)Y\ 


i<j 


FminjOi - 0j) 

{xi — axj){xj — axi) ’ 


(53) 


where x* = e®*/” and a = The functions Qkixi, ■■■,Xk) are symmetric in all variables and 

have no poles on the physical sheet. 

This ansatz, as usual in the context of the computation of form factors of local fields (see 
e.g. IMIES]), expresses the form factors in such a way as to explicitly separate the part containing 
the poles from the part which has no singularities. In addition, the explicit presence of the 
minimal form factor and the symmetry in the variables Xj automatically gives form factors 
which exhibit the correct monodromy properties in the rapidities. In the context of twist fields, 
this ansatz was used for the first time in 1691. 


4.3.1 Kinematic and bound state residue equations 


Using (46), the kinematic residue equation with the ansatz (53) can be rewritten as (A: > 0): 

Qk+ 2 {o:xo, xo, xi,..., Xk) = xlPkixo, xi,..., Xk)Qk{xi ,..., Xfc), (54) 


where Pk is the polynomial 

k 

Pfc(xo,xi,... ,Xfc) = C'fc(n) (^{xb - a^xo)(x6 - a“^xo)(x6 - /3xo)(x6 - a/3“^xo)), (55) 

6=1 


where (3 = e sn and 

Ck(n) = " a 2 (fc+i) ^ Co{n)a^^ (56) 

(k) 

Denoting a) ^ the i-th elementary symmetric polynomial on k variables xi,... ,Xfc, which can 
be defined by means of the generating function, 

k k 

x*=“Vf ^ = ]J(xi + x), (57) 

2=0 i=l 

we can rewrite Pkixo, xi,..., Xk) as 
k 

Ck{n) Y {-(y^xo)’"~°'{-a~^xo)’"~^{-af3~^xo)’"~''i-Pxo)'"~‘^a^a'>all'''ai^'^af\ (58) 

a,b,c,d=0 
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In the following we will omit the upper index when there is no confusion possible. 


Besides (54), another equation that arises from the ansatz (53) is that using the bound state 


residue equation (38). The simplest case of this equation was given in (40) and this allowed us, 


in combination with the clustering property, to fix the one-particle form factor (49). For higher 
particles, using (|4^ we find (/c > 1) 


with 


Qfc+i(xo/3 2, xo /32 , xi,..., Xk-\) = xlUk{xo, xi,..., Xk-i)Qkixo, xi,..., Xk-i), (59) 


fe-i 


[/fc(xo,xi,... ,Xfc_i) = Hk{n)Y[ixi - f3 ‘^xo){xi - 13‘^xo) 


(60) 


2=1 

k-1 




Hk{n) ^ (-/3 ^xo)^ ^ “(-/3^xo) 

a,6=0 


and 


R (r.) - . xfc _ rr 

XZA;(^) /" TT \ / • \.f/ 27 ri\( '^l(^)( 


(61) 


From the ansatz (53) it follows that Qi = 


4.3.2 Three-particle form factors 

First let us analyze the two-particle case. We have by dehnition that = Qo = {O), and 


comparing to (47), we obtain the polynomial 


Q 2 ixi,X 2 ) = {0)Co{n)a V 2 -h 


(^1 


011^2 


(O) 


((1 Q;)^fT2 - aaf) 


(62) 


It is a simple matter to verify that this is indeed in agreement with the kinematic residue 


equation (54); given Qq this is the most general solution to (54) (k = 0), as was shown in 
[6^ (in particular, the second term vanishes at xi = 0 x 2 ). Further, replacing /{O) by 

the linear combination of the zero- and one-particle form factors, Qo and Qi(xi), occurring via 
the quadratic equation (49), one can check that (62) is in agreement with (|59[). In fact, given 


arbitrary Qo and Qi(xi), the resulting expression is the unique solution to ( |M| ) {k = 0) and (59) 

{k = l). 

As was shown above, the additional condition of clustering imposes the one-particle form 
factor to take only two possible values (proportional to the vacuum expectation value), according 


to (52). For n = 1 the solution (62) is either zero (if we take the first solution in (52)) or it 
reduces to Zamolodchikov’s two particle solution for the Lee-Yang field |66] (if we take instead 


the second solution in (52)). This is in accordance with identifying the two-particle form factors 


with those of T and : F(j) ■, respectively. 

Interestingly, it turns out that the above structure subsists to higher particles: given Q 2 (xi,X 2 ) 


and Qi(xi), there is a unique solution to the kinematic and bound state residue equations (54) 
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{k = 1) and (59) {k = 2) for the polynomial Q 3 ixi,X 2 ,X 3 ). The solution has the following 
structure: 

Q3(xi, X2, xs) = Aia\a3 + A2a\al + A3aia2cr3 + A^a^ + A^a^, (63) 

where the parameters Ai are complicated functions of n but rapidity-independent. The detailed 
computation of Q 3 (xi,X 2 ,X 3 ) and the values of Ai are reported in appendix and note in 
particular that the polynomials af and (t^(T 2 have vanishing coefficients. 

Again it is interesting to consider the limit n —)• 1 of the functions Ai above. Using the 
two solutions (52), we now note that all constants vanish, Ai = 0, when we consider that 
corresponding to the operator 7~ (where ' =0 for re = 1), thus the three particle form factor 
also vanishes. On the other hand, if we consider the other solution in (52), which at re = 1 
should correspond to the field (f>, we find 


Ai — A 4 — As — 0, A 2 — —A 3 — 




I'KVfl 


231/4 


(64) 


(0) 27/2/(/^, 1)3/2’ 

and a simple computation shows that our three-particle form factor, as expected, reduces to 
Zamolodchikov’s solution |66j . 

It is tempting to use this benchmark (agreement with Zamolodchikov’s solutions) to try and 
find the general solution for higher particle numbers. However, as the three-particle case shows, 
the reduction to re = 1 occurs thanks to great simplifications. At this stage, it is unfortunately 
not obvious at all how high-particle solutions may be constructed other than by brute force 
computation. The main reason for this is the presence of two (rather than one) kinematic pole 
in the form factor ansatz (53). This leads to polynomials (xi,..., Xfc) of much higher degree 
than is the case in the standard form factor program. 


Despite the complexity of the expression (162), there are certain simplifications that can be 


used to rewrite the three-particle form factor in a form which is more suitable for numerical 
computations. It turns out that 


F3^"''(01,02,03) = f3{xi,...,Xk)ll 


Frain{di - dj 


i<j 

(Ffl')2(0)-iifi(re) 


4«sin(^)sin(|^) 


{xi — axj){xj — axi) 

3 

n^min(0i-0/), 


(65) 


i<j 


where / 3 ( 0 i, 6 * 2 , 03 ) is the function that is obtained from < 33 ( 01 , 02 , 03 ) in (162) by setting all 
terms proportional to {0)~^ to zero. In other words, when divided by Wi^j{xi — axj){xj — axi), 
all those terms simplify giving just the second summand in the formula above. This summand 
represents a kernel solution to the form factor equations, in the sense already described in 
subsection 3.2. 


Finally, note that 


lim T) 


C >|111 


a ^Ff'^( 7 o(re) ( 4 cos^ (^) X2X3 - (x2 -h X3)2) (02 - 03) 


01 


(01,02,03) = 


4sin(^)sin(|^) 

(Ffl^)2(0)-iFi(re) 

4asin(|^)sin(|^) 


(X 2 - 0 X 3 )(X 3 - 0 X 2 ) 


Fn,in(02 - 03) = 


^O| 1 ^O| 11 ^ 02 _ 03 ) 

(O) 


( 66 ) 
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where we have used the property 

Ff = {0)a-^Co{n) 




V^/ ' ' 




(67) 


which can easily be derived from (49), (56) and (61). In other words, the three-particle solution 


automatically satisfies the clustering property. This is an extremely nontrivial check of the 
validity of the three-particle solution. This situation is in contrast to that of the sinh-Gordon 
model [69], where at each particle number, the clustering property has to be imposed in order 
to uniquely hx the solution. It also follows from the result above and the cluster property of the 
two-particle form factor that 


6»i,6»2^oo \^/ 


( 68 ) 


Properties (66) and (68) are very important as they insure the convergence of the integrals (25) 
for /c = 3. 


4.4 Form factors of the fields T and ; T0 : 

In the previous subsections we have concentrated our analysis on computing the form factors of 
the fields T and ’.Tcj)'-- However, the correlators we are interested in also involve the fields T 
and : T(f> '■ thus their form factors are also required. In fact the form factors of all these helds 
are not independent from each other. We may think of T and :T4> : and of T and : T4> : as 
twist fields associated to the two opposite cyclic permutation symmetries i ^ i + 1 and i -|-1 i—i 
(z = 1,..., n, n -|- 1 = 1). From the additional symmetry under the inversion of copy numbers it 
follows that 




( 01 ,. •• ,9k) = Fl 




{Oi,--- ,0k), 


(69) 


and similarly for : F(j) : and : T(j) ■■ At the same time, as already explained in subsection 
from the non-unitarity of the theory we would expect that 


3.2 


(T) 


— 1 / 


(0i,---,0fc) = (-i)'=(r)-^F, 


— 1 ••-M/c 


{Ok, ■ 


l) 


= ,0i) (70) 


(note that (T) = (T)). These equations both define the form factors of T and impose the 
condition expressed by the last equality above on the form factors of T- We have verified that 
this is satisfied for all our solutions, and that similar equations hold for : T(t> ■■ These equations 


are the counter-part of (22) for twist fields, and show that the form factor expansion (24) is 
correct. 

Finally, another important relation which we have used in subsequent computations is the 
following identity 


F, 




{Oi,’ 


> ^k) = F^^^'"^{0i -h 27ri(/ii - 1), • • • ,0k + 2TTi{^k - 1)) < ... < Hk, 


(71) 


which allows us to express any form factor in terms of form factors involving only the particle 
living in copy 1. The same equation holds for the field :T4> ■■ 
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5 Identification of twist field operators: numerical results 


In previous sections we have provided compelling evidence for the identification of the two 
families of form factor solutions that we have obtained with the twist fields T and : Tcj) This 
evidence is based on the (highly non-trivial) fact that the one-particle and higher form factors 
of the field we identihed as T vanish at n = 1 whereas those of : T4> : reduce to the form factors 
of (j) obtained in [66]. Further evidence may be gathered by, for example, examining the short 
distance behaviour of the correlators (T(r)T’(O)) and {: T(j) : (r) : T(p '■ (0)). We must therefore 
hrst understand what the expected behaviour of such correlators should be for the theory at 
hand. 

Let us hrst consider the conformal held theory. In CFT such correlators are expected to 
converge at small distances as 


(r(r)r(0))cFT = 

(72) 

{■.T(P:{r)-.f(P-.mcFT=r-^^-^^-. 

(73) 


Indeed note that the powers above are positive for the Lee-Yang model as both c and A are 
negative (see section]^. This is of course a consequence of non-unitarity. 

In the massive theory, however, we expect that the leading short distance behaviours of these 
correlators should be described by a different power law: 

(r(r)r(O)) oc ^ 74 ) 


and 


(: Tcp : (r) : f(j) : (0)) oc 


(75) 

The reason for this is entirely analogous to the observation made in |66) regarding the 
correlator {(p{r)(p{0)). It was found that for short distance in the massive theory the leading 
behaviour of this correlator was rather than the conformal behaviour Zamolodchikov 

argued that this was due to the fact that the leading behaviour of the conformal OPE comes 
from the field (j) rather than the identity. In the massive theory the expectation value {(p) ^ 0 
and therefore the contribution to the OPE from the held (p itself becomes the dominating term 
in the short distance expansion of the two-point function. 

Similarly, it is possible to argue that the leading contribution to the OPEs of T and T and 
of : T(p : and : T(p '■ corresponds to the held (pi(p 2 ■ ■ ■ (pn where (pi represents the held (p on copy i. 
This held has dimension nA, and it is the held of smallest (most negative) conformal dimension 
that can be constructed in the n-copy Lee-Yang model. Since its expectation value is nonzero, 
it thus gives the leading contribution at short distances. Massive OPEs of twist helds will be 
discussed in more detail in section 1 

Thus, by employing a form factor expansion we may check whether the expected behaviours 
are indeed recovered from our form factor solutions. We will include up to three particle form 
factors as done in [66| . We have performed a numerical evaluation of the formula (30) including 
up to three particle form factors for the twist helds T and : T<p We conhrm with good accuracy 


that the twist helds exhibit the behaviours (74) and 


for mr <C 1. This means that (30) 


holds with xj- = Aj- — nA/2 and = A--Y(f,-. — nA/2. The tables and plots below show our 
numerical results for various n and a comparison to the exact CFT values. 
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n 

2 

3 

4 

5 

8 

10 

CFT (-4xr) 

^ = 0.3 

i = 0.756 

i = 1.15 

i = 1-52 

= 2.575 

ft = 3.26 

1-particle 

0.209643 

0.442562 

0.656773 

0.861066 

1.44896 

1.83206 

l-|-2-particles 

0.259028 

0.564549 

0.842992 

1.10754 

1.86697 

2.3611 

1-I-2-I-3-particles 

0.279487 

0.625075 

0.937636 

1.23376 

2.13554 

2.70666 


Table 1: Study of the two-point function (T(r)T(O)) of the n-copy Lee-Yang theory at short 
distances. Near the critical point we expect this correlator to exhibit a power-law behaviour of 
the form where xj- = Aj- — ^ = — ^ + • This value should be best reproduced in the 

massive theory the more form factor contributions are added. The data above show that this 
expectation is indeed met by considering up to three-particle form factors. 


n 

2 

3 

4 

5 

8 

10 

CFT {-4x,r^.,) 

A = 0 7 

i = 1.022 

— 1 35 

i = 1-68 

^ = 2.675 

1 ^ 0 04 

1-particle 

0.391185 

0.572281 

0.756341 

0.941564 

1.499823 

1.87287 

l-|-2-particles 

0.505165 

0.737822 

0.974720 

1.213628 

1.931704 

2.41539 

1-I-2-I-3-particles 

0.575841 

0.843472 

1.11533 

1.38907 

2.21646 

2.77169 


Table 2: Study of the two-point function {: T4> '■ (r) ■. T4> : (0)) of the n-copy Lee-Yang theory 
at short distances. Near the critical point we expect this correlator to exhibit a power-law 
behaviour of the form where x-T-fj,-, = — ^ = The data above show 

good agreement with CFT by considering up to three-particle form factors. 


6 Comparison with perturbed conformal field theory results 

A further consistency check of our form factor solutions may be carried out by comparing a form 
factor expansion of the correlators {T{xi)T{x2)) and {:T(t) ■ (xi) -.Tcj) '■ (X2)) to its counterpart 
in perturbed conformal field theory. 

6.1 Conformal perturbation theory and twist fields structure constants 

We may regard the action of the integrable quantum field theory as a perturbation of Lee- 
Yang CFT action by a term proportional to a coupling constant A and the CFT held (j)(x, x) of 
conformal dimension A, 

Siqft = Soft + j d^x4>{x,x), (76) 

and compute correlators by performing perturbation theory about the conformal critical point 
on the coupling A m- As is well-known, in the massive theory this coupling is related to the 
mass scale m as A oc (with a known proportionality factor [7T]). The massive correla¬ 

tors can then be obtained by using OPEs where conformal structure constants are modihed to 
structure functions of mr that can be evaluated perturbatively in A (they are convergent series 
in integer powers of A), and where vacuum expectation values, which are non-perturbative in 
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Figure 2: Graphical representation of —Axj- and for n < 11. The squares, circles and 

triangles, represent the up to one-, two- and three-particle form factor contributions. The black 
solid line represents the exact values at criticality. All curves clearly show strong linearity in n 
which is consistent with the CFT behaviour, where the coefficient of n (e.g. slope of the curves) 
approaches the CFT value as more form factor contributions are added. The agreement with 
CFT gets worse as n increases. This is also to be expected as the larger n is, the larger the 
contribution of higher particle form factors becomes (all form factor contributions are in fact 
proportional to n). 


A, are nonzero. The same type of comparison between a form factor and a perturbed CFT 
computation was carried out in [66] for the two-point function of the field (p in Lee-Yang. 

Let us now consider the OPEs of T with T and of : T 4> '■ with : Tp They involve 
only fields in the non-twisted sector (we mean by this all fields constructed by considering 
n non-interacting copies of the fields of the original theory) and by construction they must be 
invariant under cyclic permutation of the copies. Let us consider the following primary, cyclically 
invariant, homogeneous fields, composed of multilinears in the fields pi on the various copies: 
we label them by sets {ki ,..., fcj} of J different integers in [1, n] for J = 1, 2,..., re (we may 
take ki < ■ ■ ■ < kj), and take them to be 




Pkx • • • Pkj + cyclic permutations 

<Ski,...,kj 


(77) 


The symmetry factor is equal to the order of the subgroup of the cyclic replica permu¬ 

tations which preserve the sequence ki,k2, ■ ■ ■ ,kj of replica indices. That is, ^ki,...,kj is the sum, 
over all elements cr G Z„ in the cyclic replica permutation group Z^, of cr{pki • • • Pkj), divided 
by the order of the stabilizer, in Z„, of pi-i ''' Pkj- This definition guarantees that in 
every independent multilinear term, including the initial term ''' Pkj itself, appears with 
coefficient 1. The number of independent multilinears in ^ki,...,kj is n/Ski,...^kj- The symmetry 
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factors for low values of J can be written explicitly; 


•Si,k 


f 2 (n even, k = n/2 + 1) 

[ 1 (otherwise) 

f 3 (k — 1 = j — k = n + 1 — j) 

\ 1 (otherwise). 

{ 4 (k — l=j — k = p — j = n + l— p) 
2 [k — l=p — jy^j — k = n + l— p) 
1 (otherwise) 


( 78 ) 


The fields ^ki,...,kj have conformal dimensions JA. In order to have a basis of primary, 
cyclically invariant homogeneous fields of dimension JA, we need to further restrict the indices 
ki,...,kj. We may certainly fix A:i = 1, and further restrictions hold due to the residual 
equivalence relation generated by {1,..., kj} ~ {1, n+2—kj, n+l+k 2 — kj, ..., n+l+kj^i—kj}. 
More generally, in the set of all replica-index sets {ki,... ,kj}, there is a foliation by orbits, 
and a basis of fields ^ki,...,kj can be taken as fields parametrised by single representatives of 
each Z„ orbit. 

Let us give simple examples. For J = 1, we have <I>i = ^j- J = 2 the basis is 

•hi ,2 = 4 >i 4>2 + all n — 1 cyclic permutations, = (picps -|- all n — 1 cyclic permutations, etc. 
until 4*i,[n/2]+i = 4 'i 4 '[n/ 2 ]+i + all n — 1 cyclic permutations (if n is odd), or until 4*i,n/2-i-i = 
4>i4’n/2+i + all cyclic permutations up to 4>n/24'n (if n is even). In particular for n = 3, we 
have (l)i(j )2 + 4>2(t>z + only; for n = 4, we have (/>i (/>2 + (j)24>3 + and (jiiips + (j) 2 <t>A; 

etc. There is a unique field at J = n: <I>i,..,,n = 4>i4>2 ■ ■ ■ 4>n, which has dimension nA. As 
mentioned, this field is very important in non-unitary models since for A < 0 it provides the 
leading contribution (for small r) to the OPE, as it is the field of lowest conformal dimension. 

The OPEs in the massive theory can be regarded as “deformations” of the conformal OPEs 
such that the structure constants are replaced by functions of mr. Denoting by O and O any 
given pair of conjugate (i.e. whose twist actions cancel out) twist fields, it takes the form 


0{xi)d{x2) ~ r (C^^{mr)l + Cp^{mr)r‘^^<^i{x2) 

ln/2]+l 

+ C'^^''(mr)r^^4>i,fc(a:2) + ... + '’"(mr)r^”'^4>i,...,n(x2) 

k=2 

-|-Virasoro descendants, (79) 


where r := \xi — X 2 \, m is the physical mass of the Lee-Yang model. The functions 


OO 


2-2A 


(mr) = 1 + C. 


Vki,...,k) 


' (mr)^*-^ 


-2A) 




(80) 


admit an expansion in integer powers of the coupling A, hence in powers of (mr)^ and the 

constants " ’ ^ are the structure constants of the GET. In our analysis we will in fact only 
consider the leading term (the GET contribution) to these structure functions, that is, we will 
only carry out zeroth order perturbation theory whereby the mass dependence is introduced 
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through the non-vanishing expectation values of OPE fields. The analysis is still non-trivial 
because of the presence of nonzero expectation values. Note that with the definition ( |11[ ) and 
the standard definition of T and T we have the conformal normalization 




T<j}::T4>-. 


= 1 . 


(81) 


The conformal OPEs (and structure constants) of the branch point twist field T have been 
studied in several places in the literature. The most general study can be found in Appendix A 
of [71] where general formulae for the structure constants associated to the OPE of T with T in 
general (unitary) GET are given. Structure constants have also played an important role within 
the study of the entanglement of disconnected regions [za 176]. More recently the structure 
constants of other types of twist fields which arise naturally within the study of the negativity 
have been studied in [TT] [78] . However we do not know of any studies of the OPE and structure 
constants of composite fields such as : T</> :. Here we provide explicit step-by-step computations 
of the conformal structure constants and (see 

appendix [P] for details) which are proportional to one-, two-, three- and four-point functions 
of the field cj) (other structure constants would involve higher-point functions, which are harder 
to access). Our computations focus on the Lee-Yang model but could be easily generalized to 
other minimal models (for the field T the ingredients needed for such generalization are already 
provided in m)- Other structure constants and massive corrections thereof will involve higher 
point functions. The results are: 


A‘l’1 _ n _ „-4A| 

^rr - u, -n I 


1 — e" 




-4A 


for fe > 1, 


-6ATl<t> 1/1 2ni(k-l) 2^ib-l) 2mU-k) ^ ^ ^ 7 ^ 1 

= n°CT,|(l-e " )(l-e " )(1 - e - )| for j > A; > 1, 


^rr 

"tt 
'-1^1 




-6A, 

-8A 

-2A, 


27ri 27rik 27Tij 27Tip ^ „ 

{(p{e )(p[e ri j0(e ^ )(p{e ^ )) for p > j > k > 


(82) 


where k is a model-dependent function which characterizes the four-point function of fields (j) 


((/)(X1)(/)(X2)(/)(X3)(/)(X4)) = K(r?)|xi - Xi\ - X3 


I-4A 


T] = 


Xl2X2,A 

2 : 133^24 


(83) 


Other structure constants may be computed in terms of higher-point functions so that in general 
we expect 

-<i>i 


and 


with 


-2JA/i/ ^\i,/ ""‘"2 ITbbJ 

■’ =n (0(e - )0(e - )... (/)(e - )), 


^^l,k2,- -,kj _ „-2JA 


:T<t>:-.T(t>: 


= n 


K{e^ 


, e " 


k{xi,X 2 ,...) = lim \y\^^{(^{0)4>{l)(p{y)4>{xi)4>{x2) ■ ■ ■)■ 

y^oo 


(84) 

(85) 

( 86 ) 


But for the last line in (82) and for ( |85[ ) , all formulae above are particular cases of those given 
in [7i] . 
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The difficulty of calculating such terms is then reduced to the difficulty of obtaining higher- 
point functions in CFT. Such higher-point functions will also be required in order to obtain most 
massive corrections to the CFT structure constants, a problem which we will not be addressing 
in this work. 

6.2 The case n = 2 

As explained earlier, obtaining the CFT structure constants becomes a difficult problem for the 
field : T4> '■ as soon as we consider OPE terms involving products of more than two fields and 
for the field T when we consider products involving more than four fields. For this reason, the 
case re = 2 is particularly interesting as in this case the leading contribution to the OPE is given 
by the bilinear fields = 24 >i(l )2 defined earlier. The leading expansions in the massive theory 
are 

(r(r)r(O)) = + • • • (87) 


(: : (r) : : (0)) = (l + + • • • (88) 

Here the numerical coefficients arise from the total numbers of independent multilinears in <l>i 
and in d>i ^2 in the case with re = 2, which are n/Si = 2 and n/Si ^2 = 1 respectively. All 
subleading terms correspond to Virasoro descendants and massive corrections to the structure 
constants, hence are suppressed by positive powers. 

The structure constants are 


= 0 , 
rr ’ 


/=<4’1,2 _ q-SA 

^rr “ ’ 


C 






c 


^ 1,2 _ 2 


In the Lee-Yang model all the constants (|89|) can be computed. The CFT 

Cii can be found for instance in 
<p4> 


(89) 

structure constant 


cL = 

<p<p 


zr(i)ir(j)l 

5r(|)ir(|)l 


z(1.91131...). 


(90) 


The four point function of the Lee-Yang model has been studied in |79l l80l [22] . Following 
we can write the four point function as in (83) with 


'^(r/) = klM|Ti(r?)|2 + cV2(^)l'), 


(91) 


where 

Ti(r/)= ^ 2 ( 7 ?) = r?"! r?) and C = (92) 

A simple calculation then gives 

k{2) = hm \y\^^{c^{-l)m^0-)4>{y)) = -3.1802... (93) 

y^oo 
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n =2 


n =2 




n=2 n=2 




Figure 3: Zeroth order perturbed CFT versus form factor computation of the two-point function 
(: T(j) : (r) : T4> : (0)) and —(T(r)T(O)). Squares, circles and triangles represent contributions 
up to one-, two- and three-particles to the form factor expansion. For each correlator we present 
results both in linear and logarithmic scale. As expected, we see that the form factor result 
(triangles) and the CFT computation (solid line) are in relatively good agreement for small 
values of mr but quickly drift apart for larger values of mr. The range of agreement is seen 
more clearly by using a logarithmic scale, where we can directly compare the slopes of the form 
factor and CFT curves. 


Plugging these values in (89) as well as the expectation value © we obtain 

(T(r)T(O)) = rio ^1 — (4.6566...)(mr)“5^ + • • • , 


and 


(: T(j) : (r) : : (0)) = A (l - (6.2515.. .)(mr)"i -h (8.5055. ..)(mr)-t) -h • • • , 


(94) 

(95) 


which gives an approximation of the two-point function at zeroth order in perturbed CFT. 

In order to compare this with the form factor expansion, we need to hx the vacuum expec¬ 
tation values of T and : T4> Recall that we used the CFT normalization to set the coefficients 
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n 3 

of no and ra, respectively equal 1. This in principle uniquely fixes the expectation value. 
Although the resulting expectation value is not known explicitly, we may use the form factor 


expansion (32) to estimate it. From CFT the leading behaviours should be 

(4.6566. 


(r(r)r(O)) .^0 

- rvj 

(T)2 




^ — 
)r 10 




{TY 


(mr) 


and 


(: : (r) : fcP : (0)) 


4(A-A.7-^,) _ 


(8.5055...)ri ,_4 
[mr] 5. 


(96) 


(97) 


(: TcP :)2 (: Tc^ :)2 (: T<i> :)2 

We observe from Figure that the truncated form factor expansions of the two-point functions 
potentially change sign (pass by the value 0) only at short distances, at positions that become 
smaller as more particles are added. Since they approach the CFT form at short distances, 
this implies that the full two-point function never becomes zero. Since the ratios and 

(.T0.(O^T^.(o)) at large distances, they are then positive for all values of r. Hence, 

we find (T)^ < 0 and (: T(/> :)^ > 0. 


In fact, from (96) and 


we have that 


2 (4.6566...)m lo 


{TY = - 


Kt 


(: T<^ :)2 = 


(8.5055...)m "2 


K. 


-.Tcp-. 


(98) 


The constants Kj- and as expressed in (32) are necessarily positive, and the fact that two- 

point functions never become zero is related to the convergence of the series (32). A numerical 


evaluation of (32) including up to three-particle form factors yields 


Therefore 



: 1.35236, K,r^, 

1.95908. 

(99) 

{tY'' 

^ —3.443 10 

for n = 2, 

(100) 


:Y ~ 4.342m“2 

for n = 2. 

(101) 


Employing these (approximate) values in our form factor expansion we can now compare it to 
the functions (94) and (95). The results are depicted in Figures]^ 


6.3 The cases n = 3 and n = 4 

For the fields T and T it is also possible to compute the two-point function in the zeroth order 
approximation for n = 3,4. It is given by 


(r(r)'f(O)) = r 








( 102 ) 


for n = 3 (where we have used the numerical coefficients n/Si ^2 = 3 and n/5i^2,3 = l)j with 

(103) 


CW = 0, CZ]Y = 3“ 

TT ’ rr 


'*’ 1 , 2,3 _ ri-9Afj4 


= 3 

rr 
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giving 


(r(r)r(O)) = ri (^1 - (17.2221. ..)(mr)"t + (26.2893.. .)(mr)"i) +• •• for n = 3; (104) 

whereas for n = 4 we have 

(r(r)r(O)) = (^1 + + ( 4 C:^fL^ + 

^4(5^M,3^6A^^)3 ^ ... ^ (^05) 

with 


n =3 




n =4 



n =4 



Figure 4: Zeroth order perturbed CFT versus three-particle form factor computation of the two- 
point function (T(r)T(O)) for re = 3 and —(T(r)T(O)) for re = 4. Squares, circles and triangles 
represent contributions up to one-, two- and three-particles to the form factor expansion. For 
each value of re we present the same results both in linear and logarithmic scale. As expected, we 
see that the form factor result (triangles) and the CFT computation (solid line) are in relatively 
good agreement for small values of rrer but quickly drift apart for larger values of mr. 
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(106) 


^rr 

/y^l,2,3 

^Tf 


0 , 


_ A —5A 

^rr “ ’ 


r<^i,3 _ /1-6A 

^rr ’ 




‘£’ 1 , 2 , 3,4 

^Tf 


= 4-8^(0(i)</.(-l)</)(-i)0(l)), 


where we can again easily compute 


{(l){i)(l){-l)(t){-i)(l){l)) = -5.53709... 


(107) 


This gives 


{T{r)f{0)) = rx (1 - (40.7927...)(mr)“5 + (133.7569...)(mr)“5 - (120.0647...)(mr)“5 ) + 


(108) 

for n = 4. 

Like for the n = 2 case, it is possible to compare these results to a form factor expansion once 
the expectations values of T and : T4> ■ have been obtained by using (32). In the three-particle 
approximation we hnd 


giving 


and 


Kr = 2.02966 and K,r 4 >-, = 2.60713 for n = 3, (109) 

Kr = 2.89127 and K,r 4 >-, = 3.48758 for n = 4, (110) 

{Tf = 12.953 for n = 3, (111) 

{Tf = -41.5266 for n = 4. (112) 


6.4 Summary and discussion 

In this section we have studied the two-point functions {: T4> : (r) : T(j) : (0)) and (T(r)T(O)) 
by using two well-known approaches: a form factor expansion (up to 3 particles) and perturbed 
CFT (at zeroth order) calculation. Examining Figures and we can say that agreement 
between both approaches is good in terms of the range of values that the correlators take but 
not particularly good if we compare the slope and precise values the functions take at particular 
points. 

This level of agreement (and disagreement) is not entirely surprising given the expected 
range of validity of each approach: the form factors approach is eminently a large mr expansion 
and although considering contributions up to three particles should provide a relatively good 
description for small values of mr we do not expect it to be very precise for very short distances. 
Conformal perturbation theory works best near criticality, that is for very small values of mr, 
exactly where form factors should be less accurate. Besides, we have carried out perturbed CFT 
at zeroth order so the expectation is that this should really only be accurate for very small values 
of mr. Finally, a numerical comparison between CFT and form factors is only possible if the 
form factor normalization constant (that is the vacuum expectation value of the field) is known. 
In our case we can only access these expectation values approximately through yet again a form 
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factor expansion. This introduces a further error (the vacuum expectation values obtained this 
way are smaller in absolute value than their exact values) which results in an overall shift of the 
form factor points. 

Overall the results we obtain are not dissimilar to Zamolodchikov’s results |66j for the two- 
point function {(p{r)4>{0)) in Lee-Yang. Agreement with CFT was slightly better in [66] as hrst 
order corrections in perturbed CFT were also included and the exact value of was known 
from an independent thermodynamic Bethe ansatz computation. 

Despite the many limitations described above, it is still the case that agreement between 
form factor numerics and zeroth order perturbed CFT is better for some particular correlators 
than for others. We do not have a good physical explanation as to why this should be the 
case but it appears to depend on the particular functional form of the perturbed CFT curve 
obtained for each case, that is the relative weight of the various contributing terms and the 
region of values of r where the term with the lowest power of r is leading. 


7 Entanglement entropy from form factors 


7.1 Renyi entropy from form factors with mr -C 1 

In the previous section we established that although the form factor expansion is only rapidly 
convergent for mr 3> 1, it does still provide a good estimate of the short distance behaviour 
of correlators. An alternative way of testing this results is by performing a computation of the 


Renyi entropy as defined in (14). This involves also the computation of {(l)(r)(j){0)), which was 


hrst obtained in [^ (and which can be obtained from {: Tcj) : (r) : T(p : (0)) by setting n = 1). 
Figure]^ shows the results of such a computation for n = 2,4, 6 and 8. 


7.2 Bi-partite entanglement entropy of large subsystems 

In this section we will used the form factors previously obtained to study the bi-partite entan¬ 
glement entropy of the Lee-Yang model paying special attention to the region mr >> 1. 


7.2.1 Saturation 


Using (14), the entanglement entropy of non-unitary theories is 


S(r) = — lim 

n-5-i dn 


:) A{r,n) 
^^)2n B{r)'^ 


where we have used the short-hand notation 

A(r,n) := {: TP T<P : (r) : TP : {0)), 
Observe that 


B{r) := ip) {p{r)p{0)). 


lim A(r, n) = B{r) 

n—>-l 

and that lim^r-^-oo n) = limmr->-oo B{r) = 1. 


(113) 


(114) 

(115) 
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n=2 


n =4 




Figure 5: The Renyi entropy Q with ( |14[ ), subtracting the non-universal, positive, infinite 
additive contribution log-2i„ — ^ (l -T^) log(me), and evaluated in logarithmic scale using 
form factors. The form factor contributions up to one-, two- and three-particles are considered 
both for the correlators : (r) : 'T<p : (0)) and {4>{r)(j){0)). The solid line represents the CFT 

prediction log(mr) (note that for mr < 1, this is negative). All graphs show a clear 

logarithmic divergence at mr = 0 (as expected). Additional form factor contributions (also as 
expected) improve agreement with CFT. 


The expression above can be written as 


Q( \ , s , _ d A{r,n) 

Sir) = -log(me) + U — hm -— , 

^ ^ 3 ^ dn B{rY 


(116) 


The constant e is a convenient short-distance cutoff related to e by ^ log(me) = ^ log(me) -|- 
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The dimensionless, universal saturation constant is 


U = - lim 


rA'-TcP-.y [C 


'4>(j> 


n^i dn 
d 




K” 


= - lim 

n^i dn 


(117) 


where Kq was defined in 0 (equivalently (|31[)). The meaning of U is clear from a comparison 
of the small- and large-distance behaviour of the entanglement entropy: 


Sy) ~ log(me)-h t/-h o(l) 

o 


CeS 


log(r/e) o(l) 


{mr —)■ oo) 
(mr —7- 0). 


(118) 


These easily generalize to the Renyi entanglement entropy at arbitrary n, giving the saturation 
behaviour (15) with the universal saturation constant Un expressed in (16), and in particular 
U = Ui. 


7.2.2 Leading order correction to saturation 

It is easy to see that 


d A(r,n) A'(r,n) — A(r,n) log B(r) 

lim-= lim- 

n^i dn B{rY n^i B{rY 


A'[r, 1) 
B{r) 


log-B(r) 


(119) 


where A' := dA/dn. We now need to compute the objects above, all of which are given in terms 
of various two-point functions and their limits at n = 1. As we have seen in the introduction, 
both the two-point functions A(r,n), B{r) and the logarithm logil(r) admit expressions in 
terms of form factors. Here we only want to investigate the first and second order corrections 
to saturation of the entanglement so we will consider only up to the two-particle contribution 
to the form factor expansion. By doing so we have that 


A{r, n) = 1 + Ai{r, n) + A 2 {r, n) + 


where 


■f 


dd 

pr4>-.\i 

2 

^—rm cosh 6 

n 


271 

r4> :) 


TT 

(: Tcl> :) 


Ko{mr) 


is the one-particle contribution, and 


Mr,n) = ^ X] 


*J=i 

n 


/■“ /■“ deide2 

l-ooJ-oo (2vr)2 

de 




(: T4> :) 


—rm cosh Oi—rm cosh O 2 


= n 


^7-00 (2vr)2 




(: TcP :) 


Kq ( 2mr cosh - 


( 120 ) 

( 121 ) 


( 122 ) 
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is the two-particle contribution. Similarly, 

B{r) = 1 + Bi{r) + B 2 {r) ^ - 


(123) 


with 


/■“= dO 

Ff 

2 

.—rm cosh0 

1 

Ff 

Loo 2vr 

L) 

e — 

TT 

L) 


KQ^mr) = — 


3V2vr/(2«,i)2 


Ko{mr), (124) 


where the ratio was given in (51), and 
„ , , 1 roo roo 

*<’■> = 2LJ. 

/■“ cte 


Loo (27r)2 


FtiOi-e^) 


L) 


—rm cosh 9\ —rm cosh O 2 


F^ie) 


L) 


Kq ( 2mr cosh - 


(125) 


The form factor F 2 {9) was given by Zamolodchikov in [BB] and can be written as 


Ft{0) 


Trm^ i4in(6>, 1) 

8 /(ivr, 1) 


Finally 

log5(r) = Bi{r) + B 2 {r) - ]^Bi{rf‘ -\ - 

Thus, we hnd that 


, d A(r,n) 
hm -———-— 
n-5-i dn B{r)'^ 


AL, 1) - Bi{r) + A2{r, 1) - B2{r) + 


Biir)A[ir,l) + --- 


(126) 

(127) 


(128) 


where the hrst two terms will give the next-to-leading order contribution to the entanglement 
entropy (i.e. the leading correction to its saturation value) and the remaining terms give the 
next-to-next-to leading order contibution. We will now analyse this expression in more detail. 
In appendix]^ we show that A'-^{r, 1) is given by 

^'i(r,l) = Biir) + - f-L- - Ko(mr). (129) 

We also need A 2 (r, 1) which is given by 

1) = ^7v:o(2mr). (130) 

This simple result was established in [1^ for all integrable quantum field theories and even 
beyond integrability [36]. Then, the expression above simplifies to 


, d A(r,n) 
hm -———-— 
n^i dn B[rF 


o ■ \ — 7 = - I Knimr) -\ — Kn(2mr) 

/(^,1)2 VvrVS 108; ov ; g 01 ; 

d9 


3/(^,l)4 ;_oo (27r)^ 


(\Frr,m(9, 1)|^ - 1^ Kq f 2mr cosh- 


13 


33V37r/(^,l)4 


Ko{mr)^ + 


(131) 
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Thus, the von Neumann entropy of the Lee-Yang model takes the form 

2 2 / 1 \ 1 

Sir) = + 


+ 


3/( 


A POO jn , \ / I 

h:i)I (W ^ 


+ 


13 


33^^/(2p^l). 
2 

r5 


(2vr)2 
Ko{mr)‘^ + ■ ■ ■ 


Z 06 ^ 

— log(me) -I- [/ — aKo{mr) - , 

15 \/2mr 

where U is the model-dependent constant (117) and 


—2mr 2mr 

0(e-2”^''(2mr)-3/2) (132) 


2mr 


+ 


a := 


/(^,1)2 Vvrv/3 108y 


1 13 \ 

- — = 0.0769782... 



6:= a/w V - 


and 




c - 




= 0.326234... 


13 


= -0.0512159. 


(133) 


(134) 


(135) 


33\/32/(^,l)4 

In contrast to results found for unitary theories [SI 06], the results above suggest that the 
leading and next-to-leading order correction to saturation of the entropy of large blocks are 
strongly model-dependent. In particular, the leading correction is proportional to the constant 
a which clearly depends on specific features of the model under consideration (that is, the 
one-particle form factor). This term is directly related to the one-particle form factor and in 
particular to its value and the value of its derivative at n = 1. The fact that both these quantities 
are non-zero for : T4> ■ is special for this field - they would have been zero if we had used T - 
and we are tempted to conclude that this phenomenon is related to the non-unitary nature of 
the model. It would be interesting to test these results numerically for example by studying the 
spin chain model considered in [HI [20l EU 03] . 


8 Conclusions and outlook 

In this paper we have provided an in-depth study of the two-point functions of twist fields in the 
massive Lee-Yang model and their application to the computation of the bi-partite entanglement 
entropy. The main tools used for our study are branch-point twist fields and the relationship 
between their correlation functions in replica theories and the bi-partite entanglement. For 
massive unitary theories this connection was established and explored in [H], and the present 
work addresses the problem for a massive non-unitary model for the first time. Representing 
the EE using correlation functions of twist fields indeed provides the only known method so 
far for performing computations of the bi-partite entanglement in massive QET models. The 
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non-unitarity of the theory has important consequences for the computation of entanglement 
and several stark differences are found with respect to the unitary case. 

The twist field T and its conjugate T considered in |44j are not the right operators to consider 
in non-unitary models when performing entropy computations. Instead, as first proposed in [33] 
for CFT, one must consider the two-point function of suitably normalized composite fields : T(j) : 
and : T4> '■ introduced in m, defined as leading contributors to the OPE T{x)4>{0) and T{x)4>{0) 
where (j) is the lowest-dimension primary field of the model. 

In the present work we find the exact form factors of T and : : up to three particles. 

Remarkably we find that the form factor equations together with the requirement of clustering 
are sufficient to entirely fix all form factors and to provide in a natural way two families of 
solutions corresponding to the two twist fields T and : T4> ■■ We also give numerical evidence 
that the resulting correlation functions agree, at short distances, with CFT results. This is 
done by numerically evaluating truncated form factor expansions, or the logarithm thereof, of 
correlation functions at short distances, and comparing with a zeroth order perturbed CFT 
computation of the twist field two-point function, in the spirit of Zamolodchikov’s work |66j . 
The CFT computation also provides some of the first general results regarding OPEs of the 
composite twist fields : T</> : and : T(p Einally these results are used to compute the Renyi and 
von Neumann entropy, in particular in the limit of large blocks, which is well described by the 
form factor expansion. For large blocks we find that the corrections to saturation are strongly 
model-dependent for non-unitary theories. This is in contrast to the very universal form of the 
leading correction found for unitary models |44[ I46j . only depending on the particle spectrum. 

It would be interesting to compare the present results about entanglement entropy with a 
numerical evaluation in the quantum Ising model with imaginary transverse magnetic field spin 
first considered by von Gehlen [191120], whose near-critical universal region is expected to be 
described by the Lee-Yang QFT |2Tl [22] . This would provide the first strong evidence beyond 
criticality supporting the conjecture that the composite fields : T4> '■ and : T4> : are the correct 
fields for representing branch points, or conical singularities, in non-unitary models. 

The ideas in the present work and in [33] , in particular the form (14) of the twisted replica 


partition function in non-unitary QFT, lead us to speculate a relation between correlators of 
composite fields in non-unitary QFT and correlators of physical fields in its “unitary coun¬ 
terpart”. Following ideas from the field of PT-symmetric quantum mechanics [231 I24| . given 
that the non-unitary theory (i.e. with a non-Hermitian hamiltonian) considered has a real en¬ 
ergy spectrum, we may infer that there must be a similarity transformation which maps the 
Hamiltonian and correlators of the Lee-Yang model to the Hamiltonian and correlators of some 
unknown unitary theory. It is tempting to propose that the operation of taking composites with 
the lowest-dimension field </> implements, up to normalization, such a similarity transformation. 
That is, we may identify the correlators of local fields in the resulting unitary theory, denoted 
by ((O(r)O(O))), with correlators of composite fields in the non-unitary model, as 


((O(r)O'(O))) := 


(: : (r) : : (0)) 

{^{r)(p{0)) 


(136) 


It would be very interesting to test this and related ideas further. 

Several future directions of research follow naturally from this work: a more detailed study of 
the OPEs of twist fields for arbitrary n is desirable, not only to better understand the properties 
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of replica CFTs but also as building blocks for perturbed CFT computations for larger values of 
n. A systematic understanding of higher particle form factor solutions is still missing as is the 
study of the entanglement entropy of excited states and multipartite regions in massive (unitary 
or not) QFT. 
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A Definition of the field : T0 : 

Consider 

n 

■■ T(f) ■■ {y) = A lim |x - y|“ V r(?/)(/>j(x). (137) 

i=i 

The power a is fixed by requiring that the limit exist, and the normalization A is determined by 
requiring conformal normalization of the resulting field. These, as well as structure constants 
studied in Appendix may be evaluated by using standard methods of CFT [72] . Correlation 
functions with twist held insertions at yi and y 2 in the n-copy model are interpreted as correlation 
functions on a n-sheeted Riemann surface ■Mn,yi,y 2 with branch points in place of the twist helds, 
and conformal uniformization to the sphere is used. For the uniformization step, one makes use 
of the conformal map 


: A4n,?;i,y2c\ {0,CX)}, g{z) 



(138) 


with 

^ Qg ^ 1 y2-yi / - 1/1 

dz ' n{z- yi){z - 2 / 2 ) \z-y 2 ) 

In order to compute A and a, we compute the following ratio of correlation functions: 


(139) 


(: T(t> : (xi) : fct> : (^2)) 

(r(xi)r(x2)) 


= 


lim \xi — yi\ 


\x2-v2V Yl 

ii j2=i 


{T{yi)T{y2)(t>h{xi)(t)j2{x2)) 

{nyi)T{y2)) 


(140) 


The new ratio of correlators involved is interpreted as a correlator of ^j^{xi)(j)j^{x 2 ) on the 
Riemann surface Ain,yi,y 2 ^ and can be computed by using the conformal map above to relate 
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them to correlators in the complex plane. Thus 

(: Tcf) : (xi) : f (j) : (X2)) 

(r(xi)r(x2)) 

n 

= lim \xi - yi\^\x 2 - y 2 ndg{xi)\^^\dg{x 2 )\^^ V {(t){e^g{xi))(t){e^g{x 2 ))) 


2Trij2 


Vi^Xi 


ilU2 = l 


= \A\^n lim \xi — yi\ ^ ^\x 2 — y 2 \ ^ e ^ g(^xi) — e " g{x 2 ) 

Vi^^i I 

J1J2 = 1 


-4A 


= lim |xi — ^\x 2 — y 2 \°' 


Vi^^i 


X 2 - yi 


X2 - y2 


= ^^\xi — X 2 \ " lini {\xi — yi\\x 2 — y 2 \)°'^‘^^^" • 

Hence we must set 


Vi^Xi 


a = 2A (1 -) , A = n 


,2A-1 


We find the known result m that the dimension of : T(/> : is 

A 

H-7 

n 

and we have the correct normalization 

: T(j) : (xi) : T4> : (^ 2 ) ~ l|xi - X2\~^^'-'^^'-■ 


(141) 

(142) 

(143) 

(144) 


B Large n expansion of the one-particle form factor contribution 

The one particle contribution to the powers xj- and X: 7 -^: as defined in section listed in the 
tables 1 and 2 may be computed exactly as is simply given by the function (see eq. (@) 


'^\pO\l\2 


TT 


(145) 


for O = T or O =: Tcj) It is easy to show that this function admits a large n expansion in 


powers of 1/n starting with a term linear in n. Recall the expressions (52). Combining those 


with (44) we can rewrite the expectation values as 

2 


^0|i 

1 




{O) 

For large n we find that 


27rsin^sin^/(i7r,n) V VSn 


cos 




\6nJ J 


(146) 


sin 3 ^ 


27rsin^sin^V \3nJ 


^ . 2 / 2re 1 IStt^ _c, 

cos I — I + 2sm^ ( — I I = — +—h __ + 0{n ®), 


6n// TT^ 18n 9720n^ 


(147) 
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and 




27rsinj|^sin^ 




2 


2n 7 1737r^ _5 

^ ~ 9720n3 


(148) 


We now study the expansion of /(ivr, n) ^ for n large. It is possible to show that /(ivr, n) ^ = 
fkin) given by the following expression 


r (fcn + I) r {kn + I) r {kn + |) T ((fc + l)n - |) T {{k + l)ra + g) T ((fc + l)n + 
r (fcn + I) r {kn + I) r {kn + |) h ((fc + l)n - g) T ((fc + l)n - g) T ((fc + l)n + 5 ) 


It is easy to see that the leading contribution for n-large comes from the n-independent part of 
the fc = 0 term in the product. We have that 


fo{n) 

h{n) 

f2{n) 

h{n) 

h{n) 

h{n) 


and so on. This gives 


n 


VT 


(: Tcl> :) 


i r(i)r(|) 

2r(|)r(i) 

1 - ' 

1 - 
1 - 


1 - 


1 


1 


+ 


144n2 288n3 

13 19 

+ 


36n2 36n3 

+0{n-^), 

+ 0{n-*) 


+ 0{n 


-4' 


1296n2 7776n3 

25 37 _4, 

+ __, + 0{n 


5184n2 62208n3 

61 . _4. 
~ 14400n2 288000n3 ^ 

61 91 _4 

~ 32400n2 972000n3 ^ 


ir(i)r(i) /2n , 1 


—3- H— — 


197 


2r(|)r(|) n\ 24007r2 18 

(0.186944...)n - (0-0^35792...) ^ 
n 


+ - +0{n-^) 


(150) 

(151) 

(152) 

(153) 

(154) 

(155) 


(156) 


n 


TT 


pTll 2 

(T) 


i r(|)r(|) f 2 n 
2 r(|)r(i) U" 

(0.186944...)n- 


1 / 197 7 \ 

~ n V24007r2 18 J 

(0.366434...)^^ 

n 



(157) 


where the coefficient of 1/n^ is not exact but has been obtained by considering contributions up 
to /s (the next term would have given a correction of 2 x 10“^ to the coefficient). 

The expansion above indeed shows that the structure of the one-particle form factor contri¬ 
bution closely matches what is expected from CFT since 


— 4x 


_ 22 
■.T4>: - ^ 




n 1 
3 1^ 


(0.33333...)n-h 


(0.06666...) 

n 


(158) 
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and 


— 4 x 7 - 



2 n 

~5 


n 11 
3 ~ 1^ 


(0.33333...)n 


(0.73333...) 

n 


(159) 


Indeed, comparing coefficients we find that the one-particle contribution provides 56% of the 
coefficient of n for both T and : T4> 65% of the 1/n coefficient for the field : T(j) and 

49% of the same coefficient for the field T. In short, the one-particle form factor provides a 
very substantial contribution to the two-point function of twist fields, both for short and long 
distances. 


C Computation of three-particle form factors 

There are two recursive relations for the symmetric polynomial Q^: 

(53(0x0, xo,xi) = XoT’i(xo,xi)( 5 i(xi) 

= Ci{n)x^{xi — a^xo)(xi — a“^xo)(xi — /3xo)(xi — a/3“^xo), (160) 


Q?,{xf)l3 Ta^o/32,a;i) = xlU2{xo,xi)Q2{xo,xi) 


= —Hi(n)x, 


^0 


' a(F^I^)2 \ 

{0)Coin)xoXi H-— ((1 a^)xoXi - a(xo -h xj)) j 


x(xi - /3 ^xo)(xi - /3^xo). 


(161) 


The recursive equations imply that Q 3 is a symmetric polynomial of degree 6. Thus its most 

(3) 

general form in terms of elementary symmetric polynomials (Tj = cr- ' is 

(53(xi,X2,X3) = Aiafas -h + A 3 aia 2 a 3 + A^a^ + A^al + A^al + Aia\a 2 . (162) 

The constants Ai with i = 1,..., 7 are found to be: 

oFf l\(To(n) cos2 ^ - a{ 0 )-^Hi{n) cos^ ^) 


^6 — M — 0 ? Ai — A 4 — 


Ao = 


oFf I^Fi(n) - Co(n)) 


sin ^ sin i 


^sing^sinf^ 


A3 = — 




a(To(n)Ff (5 cos ^ + 4 cos ^ + 2 cos |^ + 6 cos + cos ^ + cos ^ - cos 


4 cos ^ sin sin 


a^{ 0 ) ^(Ffl ^)277i(n) (ll cos ^ + 6 cos |^ + cos |^) 


4 cos ^ sin sin 


A 5 — 


a(To(n)Ff (2 cos ^ + if (cos ^ + cos ^ + 2 cos - cos |^ + cos ^ - cos 


4cos^sinj|^sin|^ 


a^{ 0 )-\Ff^ fH 3 {n) (2 cos ^ + l)^ 


4 sin ^ sin 


(163) 
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D Conformal structure constants of twist fields 


In this appendix we present detailed computations of the conformal structure constants 


and 


:Td::T<p-. 


for different choices of the local field O. These structure constants are used in 


section in zeroth order perturbed CFT computations. The general strategy relies upon the 
fact that correlation functions of twist fields in CFT may be computed in two different ways: 
on the one hand we may treat the twist fields as standard local fields in the n-copy model on 
the manifold Ain,xi,x 2 (as defined in |138[ ); on the other hand we may conformally map the 
correlation function to the complex plane by using the map (138) thus expressing it in terms of 
correlation functions of other local fields (e.g. with no twist field insertions). Below we present 
many examples of this approach. 


D.l Structure constants involving the fields T and T 
D.1.1 The structure constant 

TT 

The CFT structure constant may be computed as follows. We may select out the term 

proportional to <hi = 4>j ™ the OPE of T{xi)T{x 2 ) by evaluating the three point function 

below, where a single field (/ii is inserted: 


(T(xi)T(x2)(/)i(x3)) 

(r(xi)r(x2)) 


Cp^\xi - X2p^ ^((/)j(x2)</>l(x3)) 

= ~ {(t>{x2)(t){xz)) 

= Ct\.\xi - X2\^^\X2 - xX\~^^ 


(164) 


On the other hand, we identify the ratio of correlators on the left-hand side as a correlator of 
4 >i{x 3 ) on the manifold M.n,xi,x 2 ^ and use the conformal map g to relate this to the one-point 
function {(j){g{x 3 )) on M^. Since this one-point function is zero in the complex plane, we have 
that, in general 


C% = 0 . 
rr 


(165) 


D.l.2 The structure constant 

rr 

The third term in the OPE T{xi)T{x 2 ) contains bilinears of the fields (/)j, with the constraint 
that they be cyclically symmetric. The only possibility are the fields 4>i ^ defined earlier, with 
k = 2,..., [n/2] -|- 1; the restriction on k is to avoid over-counting, as = ^i^n-k+ 2 - The 
coupling may be computed exactly as in the previous subsection. We consider, for 

some k G {2,..., [n/2] -|-1}, the following ratio of correlators, which we evaluate using the OPEs 
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in order to extract the structure constant: 


(T(g:i)T(x2)0l(x3)(/)fc(x4)) 

(r(xi)r(x2)) 


[t]+l ^ 

i =2 

= C'^:^'‘|xi - X2|^^(^>l,fc(x2)(/>l(x3)(/)fc(x4)) 

= - X 2 |^^(|X 2 - X 3 IIX 2 - X4|)“^^. (166) 


In the last step, we have used the fact that, by definition, every independent bilinear in $ 4 ^^ 
occurs with coefficient 1. We can then evaluate this explicitly by conformally mapping to the 
complex plane: 


(T(xi)T(x2)(/>l(x3)(/)fc(x4)) 


= l'95'(x3)P^|95((x4)|^^(((>(e‘'""5((x3))())(e‘'-'^3(x4))) 


2A/ 


(r(xi)r(x2)) 

_ n~^^|x2 - xi|^^(())(e^g(x3))(/>(e^g(a:4))) _ 

|X3 — XlP^*^^“»^|x3 — X2p^*'^’''"^|x4 — XlP^*^^“»^|x4 — X2p^*'^’''"^ 


n ^"^ 1 x 2 — xi|^^ 


27ri 

e ^ 


X 3 -X 1 

X 3 —X 2 


2 'Kik 

— e 


X4^—Xi 

X4^—X2 


-4A 


— »^|X3 — X2p^*'^’''"^|x4 — »^|X4 — X2p^*'^’''"^ 


-4Ai i4Ai — 

XI^X 2 n ^^|X 2 - Xip^|e « -e " 

rv_/ — 


I-4A 


|X3 - X2|"^^|X4 - X2|^^ 


(167) 


(where the power functions are on their principal branch), thus, comparing both formulae we 
find 


^rr 


_4A|, 2^i(fc-i) 

n 1 — e 


(168) 


D.1.3 The structure constant 

We consider now the next correction to the OPE of T and T, involving the fields ‘hi,fc,j with 
k > j > 1. Again, the ranges of k and j must be further restricted in the OPE in order 
not to overcount the fields. We do not need to discuss this in general; we just note that in 
both cases n = 3 and n = 4 there is a single field to count, 4 >i^ 2,3 = 4>i4>24>3 (for n = 3) 
and <hi, 2,3 = 4’i4>24>3 + <(' 2 </> 3<('4 + 4’34’44'i + 4>44>i4’2 (for n = 4). As usual we first consider the 
consequence of the OPE, 

(T(xi)T(x2)(/>l(x3)(/>fc(x4)((>j(x5)) 

(r(xi)r(x2)) 

X1^X2 - X2|®'^(<hl,fcj(x2)<(>l(x3)())fc(x4)(/>j(x5)) 

= C.f:^'‘’"|xi - X 2 |®'^|X 2 - X 3 |“'‘^|X 2 - X4|“^'^|X2 - Xsl”^^. (169) 
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We then perform the calculation of the correlation function by mapping to the sphere, 


X 


Xl—^X2 

r\j 


{T{xi)f{x2)(l)l{x3)(j)k{x4)(j)j{x5)) 

{T{x4)f{x2)) 

\dg{x3)\‘^^\dg{x4)\‘^'^\dg{x5)\‘^‘^{(l)ie^ g{x3))(j)(e^ g{x4))4)ie^ g{x5))) 

_I I fi A / / / / \\ I / 27rifc , , , , 2Trij 

_ra~°^|x 2 - xi|°"^((/>(e » ff(x 3 ))(/>(e n g{x 4 ,))(p{e ^ g[x^))) _ 

(|X 3 - Xi\\x4 - Xi\\x 5 - Xi|)^^(^“")(|x 3 - X2\\X4 - X 2 \\X 5 - 0 : 2 1 


) 

-2A 


( 170 ) 


C^^\x2 - xil^^n 


- Xi\\x 4 - Xill^s - "^(|X 3 - X2IIX4 - X 2 \\X 5 - ^ 21 


13^3 


n 


27Ti(k — l) 27ri(j — l) 2TTi(j — k) \ 

g{x3) - e 5(3^4)115(3^3) - e s g{x5)\\g{x4) - e s c/(x5)|j 

- x.mi - e^)(l - e^)(l - e^)|-^^ 


ka — X 2 \^^\X 4 — X2|^^|X5 — X 2 y^^ 


thus 


-6A/^<p in 2^»(fc-l) 27ri(j-l) 2niU-k) ^ ^ 

C^^’^=n CJ^|(l-e " )(l-e n )(l-e n )| 


( 171 ) 


D.1.4 


The structure constant 


^Tf 


This may be computed as before with the final results involving now a four point function of 
fields (/>: 


j',p 

^rr 


= n 


o A , , ^ , 27rifc ^ , 27ri?' ^ , 27rip ^ , 

® ((?!)(e " )(/)(e " )(/>(e " ) 4 >{^ " ))■ 


( 172 ) 


D.2 Structure constants involving the fields : T0 : and : T(j) 

Computations for the fields : T(j) ■ and : T4> ■ are very similar to those performed in the previous 
subsection, once the representation © is used. Below we provide some examples. 

D.2.1 The structure constant 

As before, we compute first 

(: T(t) : (xi) : f (j) : {x2)(t)i{x3)) a:i^3;2 I _ |2A-4A.7-^. {(t)l{x2)(t)l{x3)) 

(r(xi)r(x2)) ^ (r(xi)r(x2)) 

= - X 2 |^‘^^^“-^|X 2 - X 3 |“^^, ( 173 ) 
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and we may compute the same three point function by using the conformal map (138) together 
with the definition (11) 


(: T<t> : (xi) : f (j) : (x2)0i(x3)) 

{T{xi)T{x2)) 

n 

= lim |xi — 

Jl J2 = l 


{T{xi)T{x2)(t)j^{yi)<i)j2{y2)<i>i{xz)) 

(r(xi)r(a;2)) 


= 2 lim \xi-yi\^^^^ "'^\x 2 - y 2 ?^^^ '^^\dg{yi)\^^\dg{y 2 )\^^\dg{x^)\^^ 


Vi^Xi 

2-7rrj2 


X Y1 9{yi))4>{e ""^5'(2/2))(/>(e^"'ff(x3))) 

iij2=i 

= lim Ix2-y2|"^l'-^l|5ff(x3)p^ 

yi 


(l® ""'5(yi)-e ""£/(y2)||e 5(2/1) - e^"£/(x3)| |e ""''5(52) - e^"*5(x3)|^ 

ii j2=i 


= lim ^2-521“^^!^ "l|95(a^3)|^‘^ (15(52)1^ 15(2:3)1) 

yi 

= C^J^i-X 2 \~^ lim |a 5 (x 3 )p^| 5 (a: 3 )|“^^ 

Vi^Xi 

= - X2P^1^“"1|X3 - Xi|"^^|x3 - X2\~‘^^. 

Hence we conclude 


D.2.2 The structure constant 

Again we use 




-2A 


c 


</></>■ 


■.Tcl>::f 4 >- 


4A 

n 


(174) 

(175) 


(: T(/) : (xi) : f (j) : (x2)(/>i(x3)(/)fc(x4)) 

(r(xi)r(x2)) 

- a;2|‘^'^l^"-l((^i(x2)(/>i(x3))((/)fc(x2)0fc(x4)) 

= C^^'^,^^\xi - X2|‘^'^1^“™1 (|X2 - X3I \X2 - 2 : 41 )“'^^ ■ (176) 
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We then calculate: 

(: T4> ■ (a;i) : f (p ■ ix2)4>i{x3)4>k{x4)) 

{T{x,)T{x2)) 

= lim \xi-yi\^‘^^^~^'>\x2-y2f^^^~^^\dg{yi)\‘^'^\dg{y2)f^\dg{x3)\‘^'^\dg{x4)f‘^ 

Vi-^Xi 


Y2 giy2))(i)ie^" g{x3))4>ie"Tg{x4,))) 

iij2=i 


27rij2 


= n 


4A-2„1-4A 


^ lim \x 2 -y 2 \ \dg{x3)\^^\dg{x4)\ 

Vi^Xi 


2A 


E 27ri(?' —1) 27ri(fc —j) 

{4’{0)4'{g{y2))(t){e ^ g{x3))(j){e n g{xi))) 


i=i 


-- ■ ^ ' ' |2Am„,'„ _m2A| m-4A„ / 1 


= n 


X 1 -X 2 I “ lim \dg{x3)\'^^\dg{x4)\‘^^\g{x4)\ ( 1 - e ^ , 

V 5W3)/ 

-2A{l+i) 


-4A 


a^l — 3 ^ 21 "^^^^ "^(13:3 — Xl||x 4 — X2I) "^(|X 3 — 


X2 X4 — Xl 


k\ 1 — e" 


’inijk-l) 1' (2:4 — Xl)(X3 — X2) ^ " 
(X4 - X2)(X3 - Xi) 


xi^X2 -4A| | 4 A( 1 --)| i-4A| i-4A /i 2T^i{k-i) ■ 

~ n Xl — X 2 ^ "hX3—X 2 X 4 — X 2 K I 1 — e 


/ 27ri(fc-l) \ 

( 1 — e " j 


whence we conclude 




27ri(/c — 1 ) 


where k is a model-dependent function which characterizes the four-point function 


((/)(X 1 )(/)(X 2 )(/>(X 3 )(/>(X 4 )) = K{g)\xi - X 4 I ^^|X 2 - X 3 


I-4A 


g = 


3:123:34 

3:133:24 


(177) 

(178) 

(179) 


E Computation of A[{r,l) 


We have seen that 


n 


Ai(r,n) =- 

vr 


F, 


:r0:|l 


KQ{mr). 


(180) 


(: TcP :) 

The n-dependence of this expression is contained on the one-particle form factor, so we need to 
compute 


lim — I n 
n^i an 






(: TcP :) 


f Ff\ ^Ff d f \ 

\ {(f)) ) {(j)) n™ dn I (: r(/> :) ) ’ 


(181) 




where we have used the fact that zero real part. We now compute the derivative 

above employing the formula (52) 


lim 




i \/2 


>1 dn T(f> :) 


3i/Cf,l) 




3t/(^,l)2 dn 


• (182) 
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The derivative above can be computed from the integral representation (43) to 


lim 


dfi¥,n) 


' 27Ti 

3 

dn 


.,27ri . 


sinh I sinh | cosh ^ 

-^- ^dt 

sinh t cosh | 


Simplifying we obtain, 



(183) 


= Bi{r) + 





Ko{mr), 


(184) 


with ili(r) defined in (124). 
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